NASA TM X- 53512 


NASA TECHNICAL 
MEMORANDUM 


NASA TM X- 53512 
September 7, 1966 


(ACCESSION NUMBER) 

(THRU) 

/ 

(PAGES) 

7/WX-53 51A 

(VlASA CR OR TMX OR AD NUMBER) 

(CODE) 

32 - 

(CATEGORY) 


EFFECTS OF RECENT NASA-ARC HYPER VELOC ITY IMPACT 
RESULTS ON METEOROID FLUX AND PUNCTURE MODELS 

By Charles C. Dalton 

Aero-Astrodynamics Laboratory 

GPO PRICE $_ - 

CFSTI PRICE(S) $_ 


NASA 

George C. Marshall 
Space Flight Center, 

Huntsville, Alabama 


Hard copy (HC) 

Microfiche (MF) l 

IT 653 July 65 


TECHNICAL MEMORANDUM X-53512 


EFFECTS OF RECENT NASA-ARC*HYPERVELOC ITY IMPACT RESULTS 
ON METEOROID FLUX AND PUNCTURE MODELS 

By 

Charles C. Dalton 

George C. Marshall Space Flight Center 
Huntsville, Alabama 

ABSTRACT 


Recent NASA-ARC laboratory hypervelocity impact results are used to 
derive a formula for the just-puncturable thickness of a metallic target 
sheet. The derived formula is used for an interpretative comparison of 
photographic meteor data and the puncture data from the three Pegasus 
and two Explorer meteoroid measurement satellites. The logarithm of the 
mass -cumulative influx of meteoroids per square meter per second is found 
to be a nonlinear function of the logarithm of the least-particle mass in 
grams. A cubic polynomial for this function, with -14.52, -1.406, -0.0490, 
and -0.00074 for the values of the coefficients of respective ascending 
powers, fits the data and seems most appropriate for extrapolation over 
the mass range for cometary particles. The influx of asteroidal particles 
is estimated separately. 


Ames Research Center. 


NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 



NASA - GEORGE C. MARSHALL SPACE FLIGHT CENTER 


Technical Memorandum X-53512 


September 7, 1966 


EFFECTS OF RECENT NASA-ARC HYPERVELOCITY IMPACT RESULTS 
ON METEOROID FLUX AND PUNCTURE MODELS 

By 

Charles C. Dalton 


AEROSPACE ENVIRONMENT DIVISION 
AERO-ASTRODYNAMICS LABORATORY 



TABLE OF CONTENTS 


Page 

I. INTRODUCTION 1 

II. RESPONSE FROM COORDINATION 2 

III. EXCEPTION FOR PHOTOGRAPHIC METEORS 2 

IV. CRATER DEPTH AND PUNCTURABLE THICKNESS 3 

y. FUNCTURABILITY OF PHOTOGRAPHIC METEOR PARTICLES 5 

VI. PUNCTURE FLUX ENHANCEMENT FACTOR 5 

VII. TARGET AND PROJECTILE PARAMETERS 7 

VIII. METEOROID PUNCTURE MODEL FOR EARTH ORBIT 8 

IX. MASS-CUMULATIVE INFLUX 9 

X. ASTEROIDAL PARTICLES 14 

XI. CISLUNAR SPACE 14 

XII. INTERPLANETARY SPACE 14 

XIII. ACCURACY 16 

XIV. CONCLUSIONS 17 

APPENDIX I. Reference 11. AIAA Paper No. 66-515 23 


iii 


Symbol 

0 


LIST OF SYMBOLS 


P 

d 

*> 


P 


o 


V 


c t 


e 


t 


Pt 

P 

m 

log 


x 2 


a. 

log m 
x 

m x 

E [log m] 
v, m* 


Definition 

puncture flux; mean number per second per square meter 
of surface with 2it steradian exposure 

thickness of a jus t-puncturable sheet in centimeters 

projectile or meteoroid diameter in centimeters 

flux of meteoroids of mass equal to or greater than m; 
mean number of hits per square meter per second per 2it 
steradians without earth shielding 

thick-target crater depth in centimeters 

closing velocity in kilometers per second 

longitudinal sonic velocity in infinite target of given 
material in kilometers per second 

target ductility, relative elongation in 2- inch gauge 
length 

target specific gravity or density, gm/ cm 3 
projectile specific gravity or density, gm/cm 3 
projectile or meteoroid mass in grams 
logarithm to base ten 

angle of impact of an incident particle from the sur- 
face normal 

standard deviation of log m 

angle of impact of a puncturing particle from the sur- 
face normal 

puncturable mass at constant angle of impact x 2 

expected value or mean of log m 

antilog of E[log v] and E[log m x ] , respectively 


iv 


Symbol 
k t + log p 


+ (19/54) log m 


D 


LIST OF SYMBOLS (Continued) 

Definition 
target parameter 

target material parameter 

projectile parameter 

projectile velocity parameter 

distance from the sun in astronomical units 


v 


TECHNICAL MEMORANDUM X-53512 


EFFECTS OF RECENT NASA-ARC HYPERVELOCITY IMPACT RESULTS 
ON METEOROID FLUX AND PUNCTURE MODELS 

SUMMARY 


Recent NASA-ARC laboratory hypervelocity impact results are used 
to derive a formula for the just-puncturable thickness of a metallic 
target sheet. The derived formula is used for an interpretative com- 
parison of photographic meteor data and the puncture data from the 
three Pegasus and two Explorer meteoroid measurement satellites. The 
logarithm of the mass-cumulative influx of meteoroids per square meter 
per second is found to be a nonlinear function of the logarithm of the 
least-particle mass in grams. A cubic polynomial for this function, 
with -14.52, -1.406, -0.0490, and -0.00074 for the values of the coef- 
ficients of respective ascending powers, fits the data and seems most 
appropriate for extrapolation over the mass range for cometary particles. 
The influx of asteroidal particles is estimated separately. 


I. INTRODUCTION 


A working paper [1] prepared for consideration by the NASA Meteo- 
roid Technology Advisory Working Group* in January 1966 showed that the 
puncture data for the two Explorer and three Pegasus meteoroid measure- 
ment satellites collectively indicate that puncture flux 0 is approxi- 
mately proportional to the -1.966-power of the target equivalent thickness 
p. It was suggested that if the ratio of puncturable thickness, p, and 
particle diameter, d, can be considered proportional to the 1.04-power of 
d, then the mass-cumulative influx F^ would be proportional to the -1.34- 
power of particle mass m for the puncturing particles similarly as for 
photographic meteors. That possibility is not very convincing; e.g., 

Gault [2] had concluded previously that the ratio of the thick-target 
crater depth, p Q , and the particle diameter, d, "should never be a 
strong function of the projectile size." 

After further study and consultation, a report [3] was published 
with flux and puncture models for cis lunar and interplanetary space. 

The cislunar part was distributed as a working paper [4] again to the 
NASA-MTAWG. The models presupposed: (1) The Ames (Summers 1959 [5]) 

formulation for thick- target penetration ratio p Q /d independent of d, 


"Will be referred to as MTAWG in this report. 


(2) the Langley (see Reference 6) relation between p and thin-target 
puncturable thickness, p, (3) Hawkins and Southworth's [7,8] sporadic 
photographic meteor sample with recalculated mass values [9] and Hawkins 
and Upton's [10] value for the flux of zero-visual-magnitude meteors 
and their "visibility" statistical weighting further supplemented [11, 12], 
(4) Stuhlinger's [13] thickest- target Pegasus puncture data, and (5) for 
the Explorer XVI thickest punctured target, Hastings [14] maximum likeli- 
hood value for the slope of the logarithm of the puncture flux, 0 , as a 
function of the logarithm of the target thickness, p. All of the 
reviewers expressed a preference for the models in Reference 4 instead 
of Reference 1. 


II. RESPONSE FROM COORDINATION 


Some of the suggestions which the reviewers* offered for a further 
improvement of the models in Reference 4 are as follows: (1) The photo- 

graphic meteor tie point and slope should be determined from McCrosky and 
Posen's [15] sample which is nearly ten times as large as Hawkins' and 
Southworth's [7,8] 285-meteor sample. (2) The Pegasus puncture flux 
should be divided by 0.86 to account for the laboratory-determined count- 
ing efficiency reported by Naumann [16]. (3) The puncture data from 

Explorers XVI and XXIII should be considered collectively in the puncture 
flux model. (4) The more recent NASA -ARC hypervelocity impact data and 
analysis by Denardo and Nysmith [17] and by Fish and Summers [18] should 
be used instead of Summers' [ 5] earlier results. (5) The relation 
between thick-target penetration, p Q , and thin-target puncturable thick- 
ness, p, may depend on velocity (e.g., see Reference 19) so that p Q and 
p may scale with different powers of velocity especially for laboratory 
results with lower velocity than that for meteoroid impact. 


III. EXCEPTION FOR PHOTOGRAPHIC METEORS 


Although the results for the random sample of 285 sporadic photo- 
graphic meteors published by Hawkins and Southworth [7,8] are said to 
have been determined by a more accurate method of data reduction than 
that which was used for the larger by McCrosky and Posen [15], the latter 
sample is almost^ten times as large. Therefore, for purposes of statis- 
tical analysis, Opik (see Dalton [20]) has stated a preference for the 
larger sample. The present analysis is based on the smaller sample in 
References 7 and 8 because results were available in Reference 12 for an 
appropriately weighted statistical analysis with respect to particle 
mass with mass values re-computed in Reference 9 with Opik's [21, 1958] 
velocity dependence of meteor luminous efficiency. The mass values 
tabulated for the larger sample by McCrosky and Posen [15] are said to 


"NASA-MTAWG members and their representatives. 
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follow the practice of the Harvard Meteor Project in assuming Opik's 
[22,23] older formulation for luminous efficienty. The author’s recent 
paper [11], attached as Appendix I to this report, shows how various 
statistical results are biased independently of any constant of propor- 
tionality for luminous efficiency. Therefore, in both the present and 
previous models [3], the ordinate and slope of log mass -cumulative influx 
F> (particles per square meter per second) as a function of log cometary 
particle mass m are -13.58 and -1.34, respectively, at m = 10~ * 68 grams. 
Velocity is nominally 26.7 km/sec. 


IV. CRATER DEPTH AND PUNCTURABLE THICKNESS 


Denardo and Nysmith [17] give a tabulation and plot of data for 
impact of 2017-T4 aluminum spheres into thick targets of 2024- T4 aluminum. 
The plot shows that log p /d 19 ' 18 varies quite smoothly as a function of 
log v with slope gradually diminishing to about 2/3 at the highest veloc- 
ity data-point. A tangential extrapolation to higher velocity will be 
assumed through this data point, which was tabulated from p Q /d = 1.840, 
d = 0.159 cm (1/16 in), and v = 4.361 km/sec (22,960 ft/sec). There- 
fore, by extrapolation. 


logUPo/dXl/d) 1 / 18 ] = -0.254 + (2/3) log v (1) 


for the given projectile and target materials. 

For application of equation (1) to other target materials of the 
same density and ductility, it would seem that velocity v should be 
replaced by v/C t , as in Simmers* [5] formulation, where C t is the sonic 
velocity in the target material. Bouma and Burkitt [24] give Cj. — 6.25 
km/sec for 2024-T4. Then, by rewriting equation (1) with v replaced by 
v/(C t /6.25), 


p Q /d = 10°- 277 d 1 / 18 <v/C t ) 2 / 3 , 


( 2 ) 


for targets of the given density and ductility and the given projectile 
density. 
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The formulation with respect to the thin-target puncturable thick- 
ness p is found from equation (2) by replacing p Q with p/1.59 as in [3]; 
i.e. , 


p/d = 10°* 478 d 1 / 18 (v/ C fc ) 2 / 3 , 


( 3 ) 


for targets of the given density and ductility and the given projectile 
density. 

Fish and Summers [18] show results for puncturable thicknesses of 
various targets impacted by 0.159-centimeter-diameter balls of 2017-T4 
aluminum ranging from 0.5 to 8.5 km/sec. Stating that "although a 
threshold penetration equation is not available at this time, one can 
note several factors which affect the thickness penetrated" they showed 
a tendency for: 


p/d ~ (l/e t ) l/l8 (l/p t ) l/2 v, (4) 

where e t and p t are, respectively, the ductility and specific gravity of 
the target. Equation (4) is for the given projectile with constant d, 
etc. The velocity term in equation (4) will be ignored because it will 
be assumed that v must be replaced by v 2 / 3 as in equation (3) when the 
velocity is extrapolated to values typical of meteoroid impact. 

The values for ductility e t and density p t presupposed in equation 
(3) are those for a 2024-T4 target; i.e., [18,24], e t = 0.19 and p t = 2.77. 
Then, by equation (4) the right side of equation (3) can be multiplied by 
(0. 19/e^) 1 ' 18 for targets of different ductility. But the sonic velocity 
in equation (3) is C t ~ p^ 1 / 2 for materials [25] with the same Young's 
modulus and Poisson's ratio. Then the dependence on target density in 
equation (4) is achieved by multiplying the right side of equation (3) 
by (2.77/Pt) 5 ' 6 . Therefore, equation (3) can be replaced by 

p/d = lO 0 * 807 d l/l8 (l/e t ) l/l8 (l/p t ) 5/6 (v/C t ) 2 / 3 (5) 

for spherical projectiles of 2017-T4 aluminum with density [18] p = 2.79. 

The diameter, d, and the density, p, of the spherical projectile are 
related by 


dp l/3 = (6m/jt) l/3 , 


( 6 ) 
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where m is the projectile mass in grains. For a particle of given mass 
the right side of equation (5) is proportional to the cube root of kinetic 
energy; and for a particle of given velocity, the right side of equation 
(6) is also proportional to the cube root of kinetic energy, and is 
equated to the product of d and p 1 ' 3 . Then, if the data supporting equa- 
tion (5) had involved projectiles of different density, it should be 
expected that d in equation (5) should be replaced by d(p/2. 79) 1 / 3 ; i.e., 

p/idp 1 ' 3 ) = 10 O * 65 O(dpl/3)l/l 8 6-1/18 p-5/S (v/C fc ) 2 / 3 , (7) 


for spherical projectiles. 

By substituting from equation (6) and transposing in equation (7), 
the thickness p of any target sheet just-puncturable by any projectile 
of mass m and typically cosmic velocity is 


p = 10°* 749 mi 0 / 54 e" 1 /!® p” 5 / 6 C- 2 / 3 (v cos x 2 ) 2 / 3 , (8) 

where (as in Reference 5) the velocity v for normal incidence has been 
replaced by the normal component v cos x 2 for oblique impact. 


V. PUNCTURABILITY OF PHOTOGRAPHIC METEOR PARTICLES 


By equation (8), photographic meteor particle with mass m = lo -0 * 68 
gram, nominal closing velocity v = 26.7 km/sec, and angle of impact 
x 2 = it / 4 radians would just puncture a sheet of 2024- T3 aluminum of 
thickness p = io °* 502 = 3.18 centimeters. 


VI. PUNCTURE FLUX ENHANCEMENT FACTOR 


To establish the relation between meteoroid incident flux F> and 
puncture flux, 0 , it is necessary to consider the statistical distribu- 
tion of the mass m sufficient to puncture a given structure. Such an 
analysis was given in Reference 26 but with a different formulation for 
puncturability. Only the revisions for that analysis necessitated by 
equation (8) will be indicated here. When F> is the incident flux of 
particles with mass m k m, where m is the value of m satisfying equation 
(8) when the nominal values v and x 2 are used for v and x 2 , then puncture 

flux, 0 , is enhanced by a factor 10^ 5 ; i.e., 

0 = 10 Ps F > . (9) 
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It was found in Reference 3 that the slope p 2 of log F > as a func- 
tion of log m varies from about = -1.34 for photographic meteors to 
about half that value for the Pegasus puncture data. It should be suf- 
ficiently accurate to consider f3 5 as a constant determined from a con- 
stant p 2 = -1 and statistical independence of mass m and velocity v. 

By equation (8), and a iog v = 0.18 from Reference 26, 


log m = (54/19) [-0.749 + log(pe^ 18 p®/® C^ 3 )] - (36/19) log(v cos x 2 ) , 

( 10 ) 


a 


log m x 


0.34, 


where m x is puncturable mass at constant angle of impact x 2 . By the 
derivation in Reference 26 for E[log cos x 2 ] = -(1/2) log e, 


E[log m x ] = E[log m] - (18/19) log e - (36/19) log cos x 2 , 
m x = (m) (e" 1 ®/ 19 ) (cos x 2 )" 36 ^ 19 . 

Then, by substituting into the formulation in Reference 26 for the dif- 
ferential puncture flux, d0, corresponding to the differential impact 
angle dx 2 , 

d0 = 6.94 F^cos x) 36/ ^ 19 sin 2x dx, (11) 


where the subscript has been dropped from x in shifting attention from 
the incident particles to those that puncture. Then, by integrating 
the right side of equation (11) over the interval 0 § x ^ it/2, 


0 = 3.56 F > = 10°- 55 F . 


( 12 ) 
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The probability density function f(x) of the angle of impact, x, for the 
population of meteoroids which puncture a target exposed to an isotropic 
flux is, by equations (11) and (12), 


f(x) = (d0/dx)/0 = 1.95(cos x) 1 * 89 sin 2x. (13) 


Although the puncture flux enhancement factor in equation (12) is 
somewhat smaller than was found in Reference 26, being 3.56 instead of 
5.05, the probability density function f(x) for the angle of impact of 
puncturing meteoroids in equation (13) is not appreciably different from 
that which was found in Reference 26. Therefore, the derivation in Refer- 
ence 26 for the most probable radiant of a particle puncturing a vehicle 
of known attitude is not compromised by the new formulation for punctur- 
ability by equation (8). 


VII. TARGET AND PROJECTILE PARAMETERS 


The puncture data from meteoroid measurement satellites is for 
target sheets of different thickness and of different materials. The 
puncture flux must be determined as a function of a target parameter 
which depends not only on the thickness but also on the other signi- 
ficant parameters of the material. The logarithm of equation (8) can be 
transposed to equate target and projectile parameters by 


k + log p = k + (19/54) log m, (14) 

t P 


where the material parameter k t is 


k t = -1.360 + (1/18) log e t + (5/6) log p t + (2/3) log C t (15) 


and where the 


projectile velocity parameter 


is 


k = -0.611 + (2/3) log(v cos 
P 


x 2 ). 


(16) 
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The constant terms in equations (15) and (16) have been mutually 
adjusted so that the target parameter, k t + log p, vanishes at the 
puncture flux tie point for photographic meteors. When meteoroid mass 
and velocity are assumed to be statistically independent and either the 
flux is isotropic or the target orientation is random, then nominally, 
by equation (16), equation (14) becomes 


k + log p = 0.240 4- (19/54) log m. 


(17) 


VIII. METEOROID PUNCTURE MODEL FOR EARTH ORBIT 


In the cubic polynomial hazard model from Reference 3, 


log 0 = a Q + a 3 .(k t + log p) + a 2 (k t + log p) 2 + a 3 (k t + log p) 3 , (18) 

it is necessary to re-evaluate the constants a Q , a 3 . Because the 

target parameter, k^ + log p, vanishes nominally for photographic meteors 
with incident flux, F> = 10" 13 * 58 , and because the puncture flux enhance- 
ment factor in equation (12) is 10°* 55 , the value of a Q in equation (18) 
is 


a Q = -13.58 4- 0.55 = -13.03. (19) 

Also a x in equation (18) is the value of d(log 0)/d(k t + log p) where 
k t + log p vanishes. But for that condition, one also has from the 
photographic meteors d(log F)/d(log m) = -1.34. Then, by equations (12) 
and (17), 

a x = -1.34(54/19) = -3.81. (20) 

The remaining two constants a 2 and a 3 in equation (18) can be 
determined from the data in Reference 3 by satisfying the centroid of 
the puncture data for the thickest target on each of the three Pegasus 
satellites and the centroid of all of the puncture data for the four 
relatively thinner punctured targets on Explorers XVI and XXIII. Then, 
when the puncture flux is divided by 0.86 to account for the counting 
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efficiency reported by Naumann [16], the tie point for the Pegasus targets 
is, for the first 99 punctures, 


(log 0, k + log p) = (-7.106, -1.883). (21) 

t P 

Also, the tie point for the Explorer targets is, with the first 104 
punctures , 


(log 0, k fc + log p) £ = (-5.355, -2.633). (22) 

The values of the coefficients a 2 and a 3 which satisfy equations (18) 
through (22) are 


a 2 = -0.384 

(23) 

a 3 = —0.017. 

(24) 


Therefore, by equations (18) through (20), (23) and (24), the hazard 
model, with values compiled in Table 1 and illustrated in Figure 1, is 


log 0 = -13.03 - 3.81(k t + log p) - 0.384(k t + log p) 2 - 0.017(k t + log p) 3 . 

(25) 


The results for a 2024- T3 target sheet are also given in Table 1 and 
Figure 3. 

IX. MASS-CUMULATIVE INFLUX 

By equations (12), (19), and (25), the mass cumulative influx F> 
can be calculated from 

log F> = -13.58 - 3.81 [0.240 + (19/54) log m] - 

- 0.384[0. 240 + (19/54) log m] 2 - 0.017 [0.240 + (19/54) log m] 3 

( 26 ) 
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log F > = -14.52 - 1 .406 (log m) - 0.0490(log m) 2 - 0.00074(log m) 3 . (27) 

This result is illustrated in Figure 2 for the values compined in Table 1. 


TABLE 1 


Earth-Orbit Hazard-Model Values 


k fc + log p 

log p 

2024-T3 al 

log m 

Cometary Particles 

Asteroidal Particles 

log 0 

log F > 

log 0 

log F > 

B 

-4.098 

-13.76 

- 1.97 

- 2.52 

- 1.63 

- 2.17 


-3.898 

-13.19 

- 2.25 

- 2.80 

- 2.19 

- 2.74 

- 4.20 

-3.698 

-12.61 

- 2.54 

- 3.09 

- 2.76 

- 3.32 

- 4.00 

-3.498 

-12.05 

- 2.85 

- 3.40 

- 3.33 

- 3.88 

- 3.80 

-3.298 

-11.48 

- 3.16 

- 3.71 

- 3.90 

- 4.45 

3.60 

-3.098 

-10.91 

- 3.50 

- 4.05 

- 4.47 

- 5.02 

- 3.20 

-2.698 

- 9.78 

- 4.21 

- 4.76 

- 5.61 

- 6.15 

- 2.80 

-2.298 

- 8.64 

- 5.00 

- 5.55 

- 6.74 

- 7.29 

- 2.40 

-1.898 

- 7.50 

- 5.86 

- 6.41 

- 7.88 

- 8.43 

- 2.00 

-1.498 

- 6.37 

- 6.81 

- 7.36 

- 9.02 

- 9.56 

- 1.60 

-1.098 

- 5.23 

- 7.85 

- 8.40 

-10.15 

-10.70 

- 1.20 

- 0.698 

- 4.09 

- 8.98 

- 9.53 

-11.29 

-1-1.84 

- 0.80 

- 0.298 

- 2.96 

-10.22 

-10.77 

-12.43 

-12.97 

- 0.40 

0.102 

- 1.82 

-11.57 

-12.12 

-13.56 

-14.11 

0.00 

0.502 

- 0.68 

-13.03 

-13.58 

-14.70 

-15.25 

0.40 

0.902 

0.46 

-14.62 

-15.17 

-15.84 

-16.39 

0.80 

1.302 

1.59 

-16.33 

-16.88 

-16.97 

-17.52 

1.20 

1.702 

2.73 

-18.18 

-18.73 

-18.11 

-18.66 
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FIG, 3. MEAN CISLUNAR COMETARY METEOROID PUNCTURE FLUX 

FOR 2024-T3 ALUMINUM 



X. ASTEROIDAL PARTICLES 


As in Reference 3, the mass-cumulative influx of asteroidal particles 
is considered to be 10" 15,93 m" 1 particles per square meter per second for 
puncture-hazard purposes. This influx, with values compiled in Table 1 
and illustrated also in Figure 2 is an extrapolation of Hawkins* [27] 
model for large stones. Then, by assuming the same velocity distribution 
as for photographic meteors, it follows by equations (12) and (17) that 
the asteroidal puncture flux, 0 3 in earth-orbit can be calculated from 


log jZf= -14.70 - (54/19) (k t + log p). (28) 

Values for this result are also compiled in Table 1 and illustrated in 
Figure 1. 


XI. CISLUNAR SPACE 


The analysis of the meteoroid velocity distribution in Reference 3 
showed that the effect of the earth’s gravitational field is much less 
than had been expected. Consequently, only 25 percent of the terrestrial 
influx of cometary particles was found to be due to the gravity-effect 
factor; i.e., the gravitational enhancement factor is 4/3. According to 
Parkinson [28], the corresponding factor for asteroidal particles is 2.1. 


XII. INTERPLANETARY SPACE 


Most of the puncture hazard in cis lunar space, and in interplanetary 
space near or between the orbits of Earth and Mercury, is due to particles 
of cometary origin. This follows from Figure 1, from McCoy’s [29] pre- 
ference for Beard’s [30] estimate that the flux of cometary particles 
should be proportional to D~ 3 / 2 (where D is the distance from the sun in 
astronomical units), and from Parkinson’s [28] model for the distribution 
of asteroidal particles with respect to mass, m, and distance, D; i.e., 


> 


~ D‘ 3 / 2 . 


(29) 
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The 27 kilometers per second median lunar impact velocity, for the 
distribution developed in Reference 3, is only slightly less than the 
heliocentric circular orbital velocity for the particular heliocentric 
distance. Anyway, it would seem appropriate to assume that, at any other 
heliocentric distance, D, nominal closing velocity is proportional to 
D" 1 / 2 ; i.e. , 


v = 27D -1 / 2 , 


(30) 


for both cometary and asteroidal particles. Then, by equation (8), the 
thickness of a target nominally puncturable by a meteoroid of particular 
mass is 


p ~ D" l/3 . 


(31) 


Also, by equation (8), p 54 / 19 ~mv 36 / 19 ; i.e., together with the 
result in equation (30), it follows that the particle mass m nominally 
sufficient to puncture a particular target is 


m 


~ D ia/l9 . 


(32) 


Therefore, cometary puncture flux, 0 , is changed not only by the factor 
D -3 / 2 due to the change in the mass -cumulative flux but also by the 
factor 


D (18/19) d(log F>)/S(log m) _ D ~ 18/ 19 


These two factors combine to give 


0 ~ D-93/38. 


(33) 


These results for cometary particles, equations (29) and (33), are 
thought to be applicable between the orbits of Mercury and Jupiter, 
except that for close orbits around Jupiter it is necessary to account 
for the significant effects of the gravitational field of the planet. 
But only for D £ 1 is it safe to ignore the asteroidal particles. 
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Because of the relative sparsity of the asteroidal meteoroids in 
the terrestrial influx, pertinent information about this population is 
less substantial. Parkinson [28] preferred the following flux enhance- 
ment factors for asteroidal particles with respect to heliocentric 
distance: 1.0, 2000, 30000, and 20 for the orbits of Earth, Mars, 

asteroids (at 2.5 A.U. point of max.), and Jupiter, respectively. The 
Jovian and terrestrial gravitationally enhanced flux factors for aster- 
oidal particles were also given by Parkinson [28] as 170 and 2.1, 
respectively. 


XIII. ACCURACY 

In Reference 3 it was stated that the models which were developed 
for the influx F> and the corresponding puncture flux 0 were thought to 
be accurate to within a factor 2 probable error (normal distribution), 
or to within a factor 9 of the 3-a level of 99.86 percent one-sided con- 
fidence. The estimate of the accuracy was intended only as a quantita- 
tive expression of engineering judgment. Contrarily, the author has 
been told both that there is no scientific basis for (or value of) any 
quantitative distribution of intensity-of-belief and that the model for 
puncture flux is accurate to within a factor less than 2 even at the 
3-a level for half-mil (0.00127-centimeter) 2024-T3 aluminum foil. Cal- 
culation for this target by the model in Reference 3 gives 0 - 2.87 x 10" 4 , 
but equation (25) in the present analysis gives 0 = 1.41 x 10“ 4 , and the 
29 June 1966 model by Naumann [31] gives 0 = 7.72 x 10" 6 . Differences in 
these results are due to alternative extrapolations from essentially the 
same data. 

The author has been advised that the puncture flux 0 in equation 
(21), for the Pegasus tie point based on the first 99 punctures in the 
thickest target, should be considered accurate to within a factor 
(99 + 99 1 / 2 )/99 = 1.10 at the l-o level or 84 percent one-sided con- 
fidence. But that would involve only a conditional probability which 
would presuppose the suitability of a set of assumptions. For example, 
it would be assumed that the target had a uniform thickness over the 
entire area and of known value. The actual thickness of the aluminum 
sheet used in the capacitor puncture transducer was given by Naumann [16] 
as 406 ± 50 microns, later said to be from the specifications. Testing 
procedures ordinarily presuppose both a consumer's risk and a producer's 
risk; i.e., there is a chance that both a substandard item is accepted 
and that a good item is rejected. Anyway, a 50-micron increment in 
thickness would correspond to an increment of 0.049 in the value of the 
target parameter, k t +log p, in equation (21). And, by differentiation 
of equation (25), this would correspond to a 33 percent reduction in the 
puncture flux. 


MSFC participants deliberating on a meteoroid model for Project 
Able, August 1966. 


16 



The present analysis is thought to be more accurate than the one in 
Reference 3; but it seems prudent to retain the estimate that the influx 
F> by equation (26) and the corresponding puncture flux 0 by equation (25) 
are accurate to within a factor 2 probable error. 


XIV. CONCLUSIONS 


When the laboratory hypervelocity impact results reported by Summers 
[5], by Denardo and Nysmith [17], and by Fish and Summers [18] are inter- 
preted collectively with some restrictions for extrapolation to cosmic 
impact velocity, a relation can be anticipated between the projectile 
mass, m, in grams and the thickness, p, in centimeters of a just-punc- 
turable metallic target. The present state-of-the-art favors equating 
projectile and target parameters by 


k fc + log p = kp 4- (19/54) log m, (14) 

where the material parameter k^ is 

k fc = -1.360 + log (eJ /; 18 p| /e C| /3 ) (15) 

and where the projectile velocity parameter kp is 

kp = -0.611 + log(v cos x 2 ) 2 / 3 . (16) 

The other parameters in equations (15) and (16) are, for the target, 
ductility (relative elongation) e t , specific gravity p t , and sonic 
velocity Ct (km/sec), and for the projectile, velocity v (km/sec) and 
angle of impact x 2 from the normal. 

The statistical distributions of velocity v and angle of impact x 2 
are such that a target sheet with a 2jt-steradian exposure to the meteo- 
roid influx will have a puncture flux, 0 , exceeding the mass -cumulative 
influx, F» by a factor approximately 3.56; i.e., a particle with mass, m, 
and ordinary values of v and x 2 will just puncture the target, and F> 
presupposes m. 


17 


In earth-orbit the mean number of punctures, 0, per square meter per 
second per 2jt-s teradian exposure for cometary meteoroids can be calculated 
from 


log 0 = -13.03 - 3.81(k t + log p) - 0.384(k t + log p) 2 - 0.017(k t + log p) 3 . 

(25) 


The corresponding formula for as teroidal meteoroids is 


log 0 = -17.70 - (19/ 54) (k^ + log p). (28) 


Corresponding to equation (25) the mass -cumulative influx of cometary 
meteoroids can be calculated from 


log F> = -14.52 - 1.406 ( log m) - 0.0490(log m) 2 - 0.00074(log m) 3 ; 

(27) 


and for asteroidal meteoroids, corresponding to equation (28), 


log F^ = -15.93 - log m. 


(34) 


These results for cometary meteoroids, given numerically in Table 1 and 
illustrated in Figures 1 through 3, are thought to be accurate to within 
a factor 2 probable error, or a factor 9 at the 3a level. 

In interplanetary space the puncture flux, 0 , and mass cumulative 
flux, F>, for cometary particles are related to solar distance D by 


F >' 

^ d _3/2 

(29) 

0 ~ 

q-93/ 38. 

(33) 


but only for D ^ 1 A.U. is it safe to ignore the asteroidal particles. 
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Reference 11. AIAA Paper No. 66-515 

Velocity Dependence of Meteor Luminous Efficiency and Consequent 
Statistical Results. CHARLES C. DALTON, Scientific Assistant, Aerospace 
Environment Division, Aero-As trodynamics Laboratory, NASA-Marshall Space 
Flight Center, Huntsville, Ala. 


ABSTRACT 


A systems analysis basis of decision between Opik* s and Verniani*s 
velocity dependence of meteor luminous efficiency is given. A weighting 
function is used to transform Hawkins and Southworth f s random sample of 
285 sporadic photographic meteors into a random sample of meteoroids 
incident upon a randomly placed and randomly oriented surface at air- 
entry altitude. Each of the sample sets of mass values for the two 
luminous efficiency formulations is ranked with respect to material mass 
and partitioned at a mass value which gives subsets of comparable statis- 
tical weight. The extent to which the subsets of a set give different 
distributions of a parameter such as velocity tends to be more generally 
convincing than the numerical value of the linear correlation between 
that parameter and mass. The results favor Opik*s theory and give some 
basis, when a sufficiently large sample is used, for showing the direc- 
tional dependence of the cumulative flux with respect to a function of 
mass and velocity. 

Photographic meteor data are pertinent in the construction of models 
for meteoroid protection of personnel, equipment, vehicles, and stations 
in space because the measurable parameters have been accurately deter- 
mined for statistically significant samples. But in applying these data 
to meteoroid hazard considerations, two major difficulties are encountered, 
with the result that the analysis is less nearly direct than one might 
have expected. 

The main difficulty with meteor data in the analysis of puncturability 
is that, e.g., by extrapolation of Summers* [1] hypervelocity impact data 
to a typical meteor velocity, crater depth is a function of the density 
and kinetic energy of the particle. Thus, the luminous efficiency for 
the conversion of particle kinetic energy into meteor luminous intensity 
must be invoked. But the meteoroids not only are of uncertain chemical 
composition and physical characteristics, but also have typical closing 
velocity higher than that currently available in laboratory experiments. 
Then the luminous efficiency must be inferred from the physical theory of 
meteors . 
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Another difficulty with meteor data in the analysis of hazards in 
space is that, due to physical selectivity, a random sample of meteor 
data is produced by particles which do not constitute a random sample of 
the particle influx. Consequently, it is necessary to calculate what 
relative statistical weight the different members of the sample must be 
given so that selectivity will not bias the statistical distribution of 
the parameters of interest. 

The sample used, both in the present study and in one previously 
reported by the author [2,3] with different weighting, is Hawkins and 
Southworth's [4,5] random sample of 285 sporadic meteors doubly photo- 
graphed by the Baker Super-Schmidt cameras to a limiting photographic 
magnitude of +4 from stations at Dona Ana and Soledad, New Mexico, from 
February 1952 to July 1954. 

The main purpose in both analyses is to study the effect on the 
statistical relationships between five particle parameters - log mass, 
log velocity, arithmetic celestial latitude, heliocentric orbital eccen- 
tricity, and radiant deviation from the earth-way apex - when Opik's [6] 
and Verniani' s [7,8] formulations for luminous efficiency are inter- 
changed. The mass values tabulated by Hawkins and Southworth [5] for 
this sample are said (by Verniani [9]) to be proportional to those which 
would follow from Verniani 1 s [7,8] more recent formulation of luminous 
efficiency directly proportional to velocity. A scale translation was 
also made in the computation of mass values by Opik's [6] theory; but 
the present results are sensitive only to the velocity-dependence. 

In both the present and the previously reported study [2,3], the 
weighting function has three factors. A factor inversely proportional 
to the square of the meteor height accounts for the variation in the 
effective surveillance area. The asymmetry in the distribution with 
respect to the ecliptic plane, due to the latitude of the station, was 
reduced by weighting inversely with the fraction of the radiant circle 
of celestial latitude apparent above the horizon and by doubling the 
weight of each element with radiant north celestial latitude numerically 
greater than the numerically greatest southern celestial latitude. The 
remaining factor in the weighting function was different in the two 
analyses; being inversely proportional to the 3/ 2-power of velocity in 
that previously reported [2,3], and directly proportional to Hawkins 
and Upton's [10] inverse relative visibility g([i) in the present analysis 
where \i is the meteor magnitude above the limit of the photographic plate 
i.e. , 


g = 1 


= 10 
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= ( 0 . 71 p , - 0 . 296)" 1 


for \x > 1.85, 
for \i < 0.55, 
for 0.55 ^ p ^ 1.85. 



The previous analysis with mass values by Opik f s [6] theory (Dalton 
[2]) gave results about as one might have expected for the 285 meteor 
sample. The correlation between log mass and log velocity was only 0.01, 
the corresponding partial correlation was 0.10, and the multiple correla- 
tion for the five variables was only 0.21. But with the mass values by 
Verniani's [7,8] theory (Dalton [3]), the same correlations were -0.69, 
0.41, and 0.72, respectively. 

Meteor cosmic weight as a function of velocity v, as used by 
Whipple [11] for v > 19 km/sec, and as used for all values of v by 
McCrosky and Posen [12] and by Jacchia and Whipple [13], has a factor 
v- 2 which contains the factor v" 1 / 2 because meteor height was assumed 
to vary as v 1 ' 4 . But McCrosky and Posen [12] estimated that the uncer- 
tainty in the velocity mass law and in the number-luminosity law are 
such that corrections for these effects will probably be in error by 
at least 1 in the velocity exponent. 

There was, therefore, a reasonable suspicion that the seemingly 
anomalous correlations might represent bias due to inappropriate velocity 
weighting. But the results now to be shown for weighting as a function 
of meteor magnitude above the limit of the photographic plate do not 
support that expectation. 

Figure 1 compares Opik f s [6] meteor luminous efficiency with a hypo- 
thetical luminous efficiency proportional to v n , where the ordinate n is 
unity for Verniani f s [7,8] theory. The unit negative value for n, which 
Opik [6] suggested as an approximation, agrees very well as a typical 
value for his more detailed results, but is not assumed for the present 
analys is . 

Figure 2 shows, for each of the sets of computed mass values for the 
same meteors, the weighted correlations between the five particle param- 
eters: log mass, log velocity, radiant arithmetic celestial latitude, 

heliocentric orbital eccentricity, and radiant total elongation from the 
earth-way apex. The most striking differences are noted for the correla- 
tion r^ between the logarithms of mass and velocity (-0.83 vs -0.11), 
for the corresponding partial correlation r 12 *345 (0*35 vs 0.06), and 
for the multiple correlation r x .2345 (0*84 vs 0.20). The widely dif- 
ferent results for the correlation r 14 between log mass and eccentricity 
(-0.55 vs 0.01) seem to be consequent to the high value 0.77 for the cor- 
relation r 2 4 between log velocity and eccentricity. In the Verniani 
[7,8] theory, with luminous efficiency proportional to velocity, the 
higher velocity meteors have relatively higher luminous efficiency, and 
consequently they have a lower calculated value for particle mass. 
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Figure 3 shows that, when the sample is ordered with respect to 
the Opik [6] mass values and partitioned into subsets of equal statis- 
tical weight, the distribution of eccentricity is not very different 
for the two subsets. This result reflects the value 0.01 for the cor- 
relation r 14 between log mass and eccentricity with the Opik [6] mass 
values . 

Figure 4 shows the corresponding distributions of eccentricity in 
subsets of equal statistical weight formed by a partition of the sample 
ordered with respect to the Verniani [7,8] mass values. Although it 
would seem that these widely different distributions of eccentricity 
might be difficult to explain physically, they reflect the value 0.55 
for the correlation r 14 between log mass and eccentricity with the 
Verniani [7,8] mass values. 

Figure 5 shows, for velocity, what the preceding figure showed for 
eccentricity. These results, which reflect the value -0.83 for the 
correlation r 12 between the logarithms of mass and velocity with the 
Verniani [7,8] mass values, would seem to be even more difficult to 
explain physically. 

Figure 6, however, showing the corresponding results with the 
Opik [6] mass values, seems about what one should have expected for a 
sample of 285 meteors without any velocity weighting. These results 
reflect the value -0.11 for the correlation r 12 between the logarithms 
of mass and velocity. Therefore, it seems appropriate to consider that 
both the statistical weighting and the Opik velocity dependence of 
luminous efficiency are pertinent for a determination of the directional 
dependence of the cumulative flux with respect to a function of mass and 
velocity. 

Figure 7 shows the distribution of radiants as a function of the 
total deviation from the earth-way apex. The ordinate, the ratio of the 
probability density and the sine of the angle of displacement, is the 
relative density of radiants averaged around a circle concentric with 
the earth-way apex. The accentuated maximum in the vicinity of 90 degrees 
elongation was not expected, but it seems to be substantially established. 

Figure 8 shows the distribution of the ecliptic component of the 
elongation from the earth-way apex. The probability density functions 
for this paper, such as the one shown in this figure, were obtained by 
ordering the sample with respect to the parameter of interest and draw- 
ing a smooth curve through the vicinity of the machine-plotted point-by- 
point relative cumulative distribution. The values of the slope of that 
curve, determined at the parameter values of interest, are the values of 
the probability density function. 
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Figure 9 represents an attempt to show how, with a momentum-sensing 
impact-detection surface in earth orbit, one might infer the distribution 
of the ecliptic component of the radiant elongation from the earth -way 
apex. These results are found by deleting 30 percent of the meteors of 
lowest momentum. Figure 10 represents the corresponding results for a 
hypothetical detector with a particle energy threshold. Both results 
would seem to be inconclusive for a sample this small; but the technique 
might be useful with a larger sample. However, due to the value -0.88 of 
the correlation r 25 between log velocity and elongation, the results 
would be spurious, if there remains any appreciable bias with respect to 
velocity in the weighting function. 

Figure 11 shows the distribution of meteor radiants with respect to 
the ecliptic plane. The weighted determinations for the mean and median 
radiant deviation from the ecliptic plane are 26 and 18 degrees, respec- 
tively. 
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